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Abstract 

We compute the longitudinal hydrodynamic flow in ultrarelativistic heavy ion coUi- 
sions at ^/s = 5500 GeV by using boost non-invariant initial conditions following from 
perturbative QCD. The transfer of entropy and energy from the central region to larger 
rapidities caused by boost non-invariance is determined and the associated decrease in 
the lifetime of the system is estimated. 
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1 Introduction 



One simple scenario for treating the behaviour of QCD matter formed in the central 
region (nearly at rest in the center of mass frame) in ultrarelativistic heavy ion col- 
lisions is to neglect transverse motion and baryon number, to assume that the initial 
conditions for longitudinal motion are longitudinally boost invariant and to assume 
that the matter expands isentropically as an ideal fluid. One then obtains the Bjorken 
similarity flow [0]. The purpose of this study is to take the initial conditions essen- 
tially as given by perturbative QCD These have two characteristics. Firstly, 
they are boost non-invariant, there is no rapidity plateau but a wide gaussian-like ra- 
pidity distribution. Secondly, due to the remarkable small-x increase of the nucleon 
structure functions observed at HERA the initial energy densities are quite large, 
at LHC energies of y/s= 5.5 TeV almost 1000 GeV/fm^. This leads to rather long 
hydrodynamical evolution times (~ 100 times initial thermalisation time) and allows a 
boost non-invariant flow to develop. In the following we shall, in particular, study the 
additional longitudinal flow caused by the maxima of the initial entropy and energy 
densities at y = and the associated transfer of entropy and energy from the y = 
region to larger rapidities. 

We limit this study to LHC energies for the following two reasons. Firstly, we want 
to study longitudinal l+ld hydrodynamical effects and at LHC, due to the very large 
initial temperatures, this period lasts by far the longest. Also, during all this period 
the matter remains in the high T plasma phase. All the complications associated with 
the phase transition and the hadronic phase arise together with the need to go over to 
l-|-3d expansion Within the 1-l-ld approximation these late-time features have 

been studied in, say, 0. Secondly, at LHC energies the perturbative computation of 
the initial conditions is more reliable. At RHIC energies of v^= 200 GeV there is 
still a sizable soft component present, which makes the buildup of the initial energy 
density slower and would require separate modeling |]10|. This is even more so at the 



SPS energies of ^/s= 20 GeV, where boost non-invariant initial conditions are modelled 
by using rapidity distributions of final state particles |11| or by forcing the flow to be 
boost invariant ||12|| . 

We also emphasize that the main aim is to study the evolution of the flow under a set 
of approximations which will need corrections when applied to the physical situation. 
These include full thermalisation and validity of given initial conditions at very large 



rapidities. Also fluctuations in the initial conditions and their variation with 
transverse coordinate will have to be taken into account. 



1 



2 The equations 

For the ultrarelativistic 1+1-dimensional similarity flow it is convenient to replace 
x'^ = (t, x) by the proper time r and the space-time rapidity rj: 

r = Vt^ r7=^log^^, (1) 

2 t — X 

t = T cosh Tj, X = rsinhr;. (2) 

The variable 

i = log(r/r,) (3) 

also naturally appears. The general equation of state with one conserved quantum 
number is 

P = PiT,fi), (4) 

s = Uq-Uq = n = 3nB = -^, e = Ts-p + nn. (5) 

The aim now is to determine, for given initial conditions, the pressure p{t, x) and the 
flow v{t,x) = tanh6(t, x) from the hydrodynamic equations 

5^r^- = o, u = o,i, d,j^ = o, (6) 

where 

T'^- = (e + p)u^'u'' - pg^" , J^ = nBU^, (7) 

u'^ = (7,7f) = (cosh9,sinh6). (8) 

From eqs. (|^) and it follows that u^d^T'^" — Sfidf^J^ = Td^s^, where = su^ is 
the entropy current, so that these equations imply entropy conservation: 

d^s" = 0, 3" = su". (9) 

To express eqs. (|^) in component form it is convenient to take their components 
parallel {uyd^T^^" = 0) and orthogonal {{gau — UaUy)d^T^^^ = 0, a = 0, 1) to u^. The 
equations to be solved then become 0] 

{di + vd^)t + {t + p){vdi + d^)e = o, (lo) 

{vdi + d^)p + {e + p){di + vdr,)e = 0, (11) 

{di + vd^)nB + nsivdi + dr,)e = 0, (12) 

where 

v(r,?7) = tanh[e(r,?7) - ?7)], (13) 



2 



t = log (t/T;) 
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Figure 1: Section of space-time between two flow lines (eq. (15)) and two lines of constant 
proper time. 

with given initial conditions T = T{Ti,ri), fi = fi{Ti,ri). 

To correctly interpret the numerical results it is useful to have a concrete picture of 
the role played by the conservation laws in the flow. This is obtained by writing for 
any conserved vector 

= 1 Sxd^V = da^V^, (14) 

with dx'^ = {dt, dx) and da^ = {dx, dt) and choosing the path as shown in Fig. |l|. Here 
the horizontal lines are two lines of constant r while the vertical lines are chosen as 
flow lines r] = r7^°™(r), deflned as the solutions of 

^ = v{x{t),t) ^r^= v{t, v{r)) = tanh[e(r, r]) - r/)]. (15) 

Computing the line integral (|14D over various portions of the path in Fig. |l| gives the 
fluxes of through these portions; their total sum has to vanish. 

Take flrst = s'^. Converting the line integral in eq. (|T^) to the variables r, 1] using 
= (cosh 0, sinh 9) one finds that 

f da^s^" = I r cosh(e - r7)s [rfr/ - tanh(e - r/)rfr/r]. (16) 

Thus the entropy flux through a flow line (^) is zero {da^u^ = is an equivalent 
deflnition of a flow line) while the flow through a line r = constant is 

rm 

S{T,r]i < r] < r]2) = / dr] ts{t, rj) cosh[Q (T,r]) — r]]. (17) 

Varying now r and letting ^7i(t), ?72(r) follow flow lines, eq. (|T^ implies that the integral 
(|l^), the total entropy measured between two flow lines at a flxed proper time, is 
constant, independent of r. A similar equation holds for J'^. 
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We may also apply eq. (0) to the energy and momentum fluxes T*^^ and T^'^. Since 
they are not proportional to m'^, the result becomes more complicated. A similar 
computation gives that the flux through a flow line is: 

^ / N r psinh9 , , 

FEin <T <T2) = - dr -— — — -, (18) 

^ ^ Jri cosh[e - ?7fl°™(r)] ^ ^ 

/ N P J p cosh 6 

Fpin < r < Ta) = - / dr -f- — --.r, (19) 

^' Jr^ cosh[e -r7fl°^(r)]' ^ ' 

where the arguments of Q,p are r, r7^°™(r). Further, the energy and momentum fluxes 
through a segment r = constant are: 

E{T,r]i < 7] < r]2) = / dr/ r[e cosh 9 cosh(9 — 77) + psinh 9 sinh(9 — ■?])], (20) 

Jr]i 

rv2 

P{t,7]i < T] < 7]2) = I dr/ r[e sinh 9 cosh(9 — 77) + pcosh 9 sinh(9 — r])], (21) 

where the arguments of e,p, 9 are r, r/. Referring to Fig. |l|, the total sum of two 
contributions of type ([T8|) and of two contributions of type ( pOf ) has to vanish for T^^\ 
similarly for T^^. This implies that the energy (eq. (|20D ) or momentum (eq. (^)) 
between two flow lines is not constant but changes due to energy or momentum flow 
across the flow line: work done against expansion. 

Note the different flow-dependent factors in (|1^ and (pO]): the total energy E con- 
tains an additional boost factor cosh 9 not present for S. 

To have a still simpler view of the conservation laws, assume a similarity flow, 
9(r, r/) = Tj 01 v(t,x) = x/t and an equation of state p = p(e). Then v = and 
from eq. (|T5|) the flow lines are 77 = constant and we take them to be ±r/o. For this 
very special flow the equations of motion (|^) become 

rdre + e + p = 0, (22) 

= 0, (23) 

rdrU + n = 0, (24) 

rdrS + s = 0, (25) 

(the last follows from the flrst and third and e + p = Ts + /in) and further imply that 

T' T' 

P = p{t), e = e(r), s = s{t, t]) = -s{t„ t]), n = n(r, 1]) = -n{Ti, 1]), (26) 

r r 

i.e. p, e depend on r only while s, n can depend also on r] in such a way that Ts + fj,n 
depends only on r. The conservation laws of s^,T^^ then simplify to 

m 

dr] [Tis{Ti,r]) - T2s{T2,ri)] = 0, (27) 

/VO /•T2 
dr] cosh T] [Tie{Ti) - T2e{T2)] = 2/ drp(r) sinhr/o; (28) 
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Figure 2: The initial condition for e{Ti,ri) for ^/s = 5500 GeV. The dashed curve is the 
gaussian fit, eq. (|3l|), with a = 3.8. 

similar ones hold for J'^, T^^. Eq.(^) is clearly an integrated form of eq. (pS]) . In eq.(p8D 
the rj dependent parts factor and match on two sides of the equation. The r dependent 
part explicitly shows how the change in e(r)r is related to pdV ~ pdr. 

A further step of simplification would be to assume a massless equation of state, 
p = e/3 = aT^ + bfi'^T'^ + cfi^, a,b,c = constants. Then a solution with a similarity 
flow would be 

/r\i/3 /r\i/3 /r\^/3 

T{T,r^)=[^fj T{n,r^), ^^{r,r^) = [^fj M^.,^), = (7) (29) 

where T{Ti,ri), fi[Ti,ri) are so constrained that p{Ti) is independent of rj. This is the 
standard Bjorken flow [0] generalised by the inclusion of baryon number. Putting /i = 
finally gives the Bjorken flow. 

3 Initial conditions 

One expects the initial net baryon number to very small near y = 0; in |Q the net 
baryon number-to-entropy ratio, B/S was computed to be about 1/5000. Hence, we 
shall put the chemical potential fi = 0. Further, we are interested in the evolution of 
the system in the plasma phase, T = Tc...7Tc and choose the equation of state to be 
p = p{T) = e/3 = aT^, a = constant, s = p'(T) = {e + p)/T. According to lattice data 
for pure SU(3) [|1^] the validity of p = e/3 improves with increasing T so that the error 
is 20% at T = 2Tc and 10% at T = 3Tc. 

The initial conditions have to specify e(rj, 77) and 0(rj, 77). For the initial flow we shall 
simply take a similarity flow, Q{Ti,ri) = t], in order to study how boost non-invariance 
affects it. The initial conditions for the energy density are computed by extending the 
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calculations of to all rapidities. The initial energy per unit rapidity in the local 

rest frame equals the transverse energy of produced minijets computed from 



N dE 1 ^ , r , da^^ 1 , , 

V) = -i-TT ~ ^AA(b = 0) / dpTPT-: — -f- ' tt, (30) 
arj V Jpo dpTdy v 

where the nuclear overlap function is Taa(O) ~ A'^/{7tR\) ^ 32/mb for Pb+Pb, the 
volume per unit rapidity is V = Ti7rR\, 1/Ti = po = 2 GeV and the inclusive gluon jet 
production cross section in NN collisions is computed in 0. The results for = 0.1 
fm are shown in fig. |^ for LHC energies, ^/i= 5500 GeV, together with a gaussian fit 
to the central region 

e(r,,r/) =eoexp[-r/V(2a2)], (31) 

with a = 3.8. The distribution is seen to be quite broad. In this sense boost non- 
invariance is quite mild in the central region, even though evidently the initial condition 
has to drop faster than a gaussian at the ends of phase space. As a side remark, at 
RHIC energy \/s= 200 GeV the width parameter is a = 2.05. 

In the initial values of in fig. |^ only gluons with pt > 2 GeV are included. The value 
2 GeV corresponds to the saturation limit at LHC energies and the resulting can 
be expected to be a good estimate for the total energy density. At RHIC energies the 
saturation limit is lower, only about 1 GeV. This is so low a scale that perturbation 
theory becomes unreliable. One could also try to keep the minimum pt of partons 
included in the computation at the fixed value of 2 GeV; at RHIC energies one then 
should also include a sizable (about 50%) soft component if one wants to model the 
entire event. Modeling the soft component is certainly possible but we wish to avoid 
this phenomenological analysis [l^ in the present work. 



According to eq. ( [20| ) the total initial energy carried by the flow in the interval —r]Q < 
Tj < TjQ is, estimating the transverse area to be vri?^, 

Etotin) = nRl / ndvcoshf] e(r„ r]) ^ / d'p E—— ■ Taa{0)- (32) 

J~-rio •'PT>Po,\y\<VO " P 

Here the first part is a general relation for a flow = 77; the second part shows where 
it is computed from. Due to the boost factor cosh 77, which in the jet computation 
corresponds to the weight factor E, most of the energy is carried by flow at large ra- 
pidities. When TjQ approaches the beam rapidity ?/beam, the error in the jet computation 
grows for two reasons: firstly, baryon number will be important at large rapidities and 
secondly, physics there is not that of independent 2— 2 collisions. However, this error 
is important only at very large values of rj, where the energy density is already small. 
It thus will not affect the conclusions of this hydrodynamical study. 
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4 Analytic approximations 

We are studying modifications to the Bjorken flow and it is appropriate to ask whether 
they in some hmit can be described analytically. Consider the limit cr ^ 1, a broad 
gaussian Writing (t = log(r/rj)) 

4 - r7^ 

e{T,v) = eoexp{--t - — + ^(r)}, (33) 
e(r,r/) = r^[l + /(r)], (34) 



linearising eqs. ([TT| - |T2D under the assumptions f^g^rf/a"^ ^ 1 and integrating, one 
obtains 

f{r) = -^[l-exp(-|£)], (35) 

9ir) = -^{t-^[l-exp(-^t)]}. (36) 
The linearisation dj:f ^ demands that 

£«31og(^a). (37) 



Furthermore, for the flow lines one obtains 

log[r^«°-(r)^o] = ^{t - ^[1 - exp(-^£)]} (38) 

These equations contain the following natural expectations: 

• The energy density at rj ^ decreases faster than in the Bjorken flow due to 
energy moving to larger rj. For small times 

e(r) = eoexp[-^(t + ^ + ...)]; (39) 

• With increasing time the flow is accelerated relative to the similarity flow. For 
small times 

e(r,r/)=r/(l + ^ + ...); (40) 

• The flow lines bend outwards. For small times 

r^fl°-(r)=r^o(l + ^ + ...)• (41) 
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Figure 3: The energy density scaled by the Bjorken flow, e(r, •q)/ei*- (j jTi)'^!'^ ^ for i/s = 5500 
GeV. From top the curves correspond to r/r^ = 1, lO^/^^ 10, 10^/2, 100. 




Figure 4: Curves of constant e(T, rf) for ^fs = 5500 GeV. The curves from below at r/ = 
correspond to e = 900,117,15,2 GeV/fm^. The vertical lines are flow lines {r\ = constant 
for Bjorken flow). The thick flow lines are used to compute fig. ^. The dotted lines are the 
leftmost characteristic curve and the rightmost C~, between which the computation has 
to remain. 

5 Numerical results 



The method of characteristics is particularly suited for 1+1 dimensional hydrody- 
namical problems. First the differential equation 



— ^ = tanh[6(t, rf) 
dt 



(42) 
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where Us = tanhc^ = tanh(l/v^), defines the two famihes of characteristics. This 
is hke the equation ( [T5| ) defining the fiowhnes, but modified by the sound rapidity ±ys. 
Along each family of curves the changes of and e are related by 

de ± = 0, along (43) 
e + p 

A code integrating e, 6, starting from the initial values e(rj, i]), 6(rj, r^) = 77 by deter- 
mining the characteristic directions, and then the new values of e, by stepping in the 
characteristic directions, is easily constructed. 

Results of numerical integration of the equations are shown in Figs. for LHC 
{^/s= 5500 GeV). The energy density, scaled by the Bjorken fiow, is shown in fig. |^ as 
a function of r] at different times. Curves of constant energy density on the r, r] plane 
are plotted in fig. ^. The energy density in fig. ^ at 77 = is seen to decrease somewhat 
faster than in the Bjorken case. This is due to leakage of energy to larger rapidities, 
seen as an increase in e at large r] and described analytically in section 4. During the 
whole duration of the plasma phase, e>2 GeV/im^, or from to about lOOxj, this 
additional decrease is about 12% at LHC. 

Converted to the lifetime of the system in the plasma phase (the time it takes 
to decrease from to ec) the above result implies that the density gradient in the 
longitudinal direction decreases the lifetime by about 9%. This is opposite to the effects 
of dissipation: in the case of a rapidity plateau, the fastest decrease of energy density 
is obtained in the case of full thermalisation; dissipative effects increase the lifetime of 



the system |T0|. The longest lifetime is obtained for free-streaming expansion. 



Fig. ^ also shows flow lines, which are constant in rj in the Bjorken case but bend 
outwards in the present case. We do not show the numerical grid of characteristic 
curves, but the leftmost characteristic curve and the rightmost one C~ are plotted 
in this flgure. The numerical computation cannot be extended outside them. 

The flow 0(r, r]), scaled by the Bjorken flow rj, is shown in flg. ^ The initial flow 
Q = 1] is rapidly accelerated at large \r]\, due to increasing dr/p/p, (=rj/(T'^ for the 
gaussian parametrization) but this is already in the domain where the details of the 
model need not be correct. In the relevant region near 77 = the effects are seen to be 
small, a few %. 

As discussed earlier, the total entropy between two flow lines, eq. (p!7|) , is constant. 
However, the flow rapidity G(r, 77) changes along the flow lines changing also the ra- 
pidity interval between the flow lines. As a result the entropy of given fluid element is 
shifted in rapidity and its amount per unit flow rapidity is changed. With our initial 
conditions the shift is always to larger rapidities — higher pressure in the central region 
accelerates the fluid towards the ends — and the rapidity intervals between the flow 
lines increase. The net result is that the entropy per unit flow rapidity, which can be 
expressed as 



dS(T,r]) -r^l , \ 1 ^/-^/ \ 

— = TrR^Ts{T, T]) cosh[e(r, ^) - 



(44) 

fixed T 
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Figure 5: The flow rapidity scaled by rj, Q{T,rj)/r], for ^/s = 5500 GeV. The curves at rj = 
correspond from below to values of r as in fig. |^. 
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Figure 6: The total entropy per unit Q between the two thick flow lines in fig. |^ for ^/s 
5500 GeV. The dotted curve corresponds to r = Tj and the solid one to r = lOOrj. 



will decrease in the central region and the overall distribution will get wider. 

The result of our computation for dS/dQ is shown in Fig. |^. In t] the computation 
is extended between the two (arbitrarily chosen) thick flow lines in Fig. ^. The dotted 
curve shows the initial distribution at = 0.1 fm and the solid curve the final distri- 
bution at r = IOOtj. Note that the end points move outwards due to a combination of 
two effects: the flow lines bend outwards (Fig. ^ and the flow is accelerated relative 
to 7] (Fig. D. 

To obtain a measurable distribution we should be able to treat the hadronization 
and to fold the thermal motion of final particles with the collective motion of the fiow 
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(after expressing the entropy in terms of particle number densities). At present, we 
do not have a rehable way to estimate the effects of hadronization but if it has any 
effects on flow we would expect them to further widen the rapidity distribution as the 
hadronization proceeds from the lower density fragmentation regions to the central 
region. Thermal folding involves flow velocities effectively over two-to-three rapidity 
units. It will lead to a somewhat wider overall rapidity distribution but in the smooth 
central region the effect of folding is small. We conclude that our result of Fig. ^ gives 
the minimum change from the initial state at to the rapidity distribution of final 
particles. 

6 Conclusions 

In the case of a fully developed rapidity plateau, there is a very simple hydrodynamical 
scaling solution for the longitudinal expansion of QCD matter produced in ultrarela- 
tivistic heavy ion collisions [0]. In reality, there need not be a rapidity plateau, and 
we have studied the longitudinal flow using a rapidity distribution of the initial energy 
density obtained from computations in perturbative QCD. Then the rapidity distri- 
bution is approximately gaussian but very broad, and thus the deviations from the 
Bjorken flow are not very large. To analyze the general features of the flow we also 
derived equations for the total entropy and energy between two flow lines and gave 
approximate solutions in the limit of broad gaussians for the central region. 

The finite width of the rapidity region leads to a transfer of energy from the central 
region to larger rapidities. As a consequence the rapidity distribution gets wider and 
the energy density in the central region decays faster than for a boost invariant flow. 
Correspondingly, the lifetime of the system (the time it takes to decrease from ei to 
ec) decreases by about 10%. A similar decrease would be caused by transverse density 
gradients. On the other hand, dissipative effects would increase the lifetime. 

The longitudinal density gradient also leads to some but very small acceleration of 
the longitudinal flow: the effect in the central region is rather on the 1% than 10% 
level. In the sense our results confirm the assumptions in §], called the frozen motion 
model in |T^, that in estimating thermal effects in the central region it is reasonable 
to assume that at LHC energies the longitudinal flow velocity scales even when the 
density distributions are not boost invariant. 
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